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We show that a superconducting circuit containing two loops, when treated with 
Macroscopic Quantum Coherence (MQC) theory, constitutes a complete two-bit 
quantum computer. The manipulation of the system is easily implemented with 
alternating magnetic fields. A universal set of quantum gates is deemed available 
by means of all unitary single bit operations and a controlled-not (cnot) sequence. 
We use multi-dimensional MQC theory and time-dependent first order perturbation 
theory to analyze the model. Our calculations show that a two qubit arrangement, 
each having a diameter of 200nm, operating in the flux regime can be operated with 
a static magnetic field of ~ O.IT, and an alternating dynamic magnetic field of am- 
plitude ~ 1 Gauss and frequency ~ lOHz. The operational time Top is estimated to 
be ~ 10ns. 
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I. INTRODUCTION 

In the last few decades there has been broad interest in hope to design and 
construct a practical quantum computer. The, so to speak, machines will enable 
us to reach a domain of knowledge that was, up to now, considered unreachable or 
beyond human capabilities. They will allow us to perform calculations of enormous 
amount in a rapid and effective manner. Factoring large numbers, teleporting large 
amounts of information, and simulation of the real world, will enter our immediate 
reach and are destined to dramatically changes our lives. 

The core role of quantum computers are played by quantum bits, qubits. Many 
ideas have been raised for practical qubit realization: cavity quantum electrody- 
namics, ion traps and nuclear spins. Also, superconducting circuits of Josephson 
junctions have been proposed and experimental research is being pursued. In this 
paper we will show a possibility to implement a quantum computation scheme on a 
coupled flux qubit system. In the rest of this introduction we will give a brief descrip- 
tion of the idea of quantum computation, and also, an outline of the macroscopic 
quantum coherence theory (MQC), with which we treat our system. 

A. Quantum computation 

Quantum computers have attracted a lot of attention and an abundance of work 
has been established. But here, we will concentrate on only the main aspects of 
quantum computers. Least requirements for a quantum computer are that the 
system must retain certain properties: 

1. Ability to represent quantum information, meaning that a quantum bit must 
be able to represent not only two classical values |0) and but also a super- 
position of these two states, namely a|0) where a and f3 are c-numbers. 

2. Have a universal family of unitary transformations. Universal, here represents 
that the Hamiltonian of the system is capable of controlling the system state 
arbitrarily. In other words one must be able to reversibly transform any given 
state into another state of choice. The reversibility is a quantum mechanical 
requirement. DiVincenzo 0] showed that a controUed-not (cnot) and single 
qubit gates were universal for any n-qubit system. 

3. Have a preparable initial state. This is an obvious demand, since the initial 
state of the system is the input. However there is no necessity for capability 
to provide an arbitrary initial state, since the above requirement suffices to 
transform a particular initial state into the desired input. 

4. Have a means of measurement. One must be able to measure the probability 
amplitude of the final state, i.e. the output. Measurement of the system state 
is often the result of the qubit's couphng to a classical system, commonly 
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the environment. This process is equivalent to a projection of the superposed 
quantum state of the qubits. Although each measurement outcome is generi- 
cally random, by controlling ensembles of quantum computers one is capable 
of determining the output state. In fact, a measurement does destroy the 
superposition and acts randomly, so one must be sensitive to unwanted mea- 
surement which is a cause of decoherence. 

All quantum computers must satisfy the above four requirements, plus many oth- 
ers for optimal and efficient computing. A qubit, due to its quantum nature, suffers 
from decoherence, where the coherent superposition of quantum states is destroyed 
by noise. This is the most influential obstacle in creating a practical quantum com- 
puter. Extracting all noise sources from the system is definitely impossible, but if 
we could perform our operation before the system loses its coherence then we could 
obtain our output with relatively high fidelity. Provided, a strong requirement for 
an efficient quantum computer would be for the 'quality factor', r^/xop > 10"^ jif. 
Top is the time necessary for a single operation, and is the decoherence time: the 
time length that the system can maintain its coherence. 

The coupling of qubits play an important role in quantum computation. To 
create large scale accurate quantum computers, there must exist an efficient method 
of coupling selected qubits, and to apply transformations upon the coupled qubits. 
Although coupling within the system is permissible, often the system couples to the 
environment thus resulting in considerable noise. This is an important point. 



B. The quantum mechanics of flux qubits 

In our research, we have chosen superconducting circuits as our qubit arrange- 
ment. The advantages are that: they are relatively easy to fabricate, they can be 
measured easily, and large arrays can be effectively implemented. Many of these 
characteristics are results due to the fact of the flux qubit being a macroscopic 
device. 

However, because of its macroscopic nature an obvious question rises: Does it 
behave purely quantum mechanically? In other words, can the device be put in a su- 
perposition of two distinct states? The answer to this question is provided by Macro- 
scopic Quantum Tunneling or Macroscopic Quantum Coherence (MQT/MQC) the- 
ory. This field of study concerning superconducting loops was first developed theo- 
retically by Ivanchenko et. al. j^, and Caldeira and Leggett 0. Experiments that 
followed Bi have confirmed " fa. that the magnitude of magnetic fiux piercing 
a superconducting loop, when seen as a canonical variable, does obey the rules of 
quantum mechanics. 

We distinguish each state of the system by the direction of electric current, so 
each state is macroscopically distinct. Here, we step out of fundamental quantum 
mechanics and see that a macroscopic matter can take a state not definitely \a) 
nor \b), but a\a) + (3\b). This occurs paradoxical to the sane mind, and has been 
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a discipline of long discussion 0: can the cat be dead and alive? Although it is 
a fundamental aspect, we will not indulge ourselves with a discussion of the EPR 
paradox here. 

FlTthermore, theory 0, S 13 and experiments Q ^ have shown that tran- 

sition between these macroscopically distinct states can be induced by photons or 
dynamically alternating magnetic fields. 

Consequently, the tools of the trade for quantum computers can be provided 
with superconducting circuits. However there still has not been found proof to 
whether MQT/MQC theory can be applied to systems with more than one degree of 
freedom, each being a macroscopic value. Experiments j^, 14| have shown evidence 
that MQT does occur in the thermal regime, but also indicate that the interaction 
between degrees of freedom, resulting in a suppression of escape rates, cannot be 
ignored. In a recent experiment 0, the quantum regime MQT theory seems to 
agree well with the behavior of a system with more than one degree of freedom. 



II. DOUBLE QUBIT SYSTEM 



A. Description of model 




FIG. 1: The double qubit model. The crosses represent Josephson junctions, and the lines 
superconducting nodes. Each loop represents a qubit, and tunneling current circulates the 
circuit. 

In this paper, we have taken a superconducting circuit as is shown in Fig. d The 
model consists of six Josephson junctions and five superconducting nodes. The two 
loops of the circuit correspond to the two qubits. 



B. Static properties 

We assume that the temperature can be lowered well below Ag/^b, where 
is the gap energy of the superconductor and ks the Boltzmann constant. Allowing 
this assumption we proceed with our analysis on the basis that each node is in a 
coherent state, i.e. maintaining a single value order parameter throughout the node. 

We apply a uniform static magnetic field perpendicular to the plane the circuit 
rests in. By taking the phase difference of each junction 7j as our macroscopic 
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variable, from Josephson's Law, the tunnel current flowing through each junction 
becomes, 

4 = 4*=sin7fc, (1) 

where, is the critical current of the k-th junction. From the single valuedness of 
the order parameter, we see that the 7s must always satisfy the conditions, 

71 + 72-73 = 27r/, (2a) 
-72 + 74 + 75 - 76 = 27r/, (2b) 

where / represents the magnetic flux through the loops, measured in units of flux 
quanta ^0, and is often referred to as the frustration index. The junction energy of 
this system can be seen to be 

W(7) = UM+U^aM 

6 

= {1 - COS -fi) +Umag{l), (3) 

i=l 

where Ei represents the junction of energy of the i-th junction, and Umagij) is the 
magnetic energy of the system. This is an analogous extension of the washboard 
model to multi-dimensions. According to MQT/MQC theory, the whole system can 
be regarded as a free particle with position 7 moving within a potential U. 

For the kinetic energy of the system we take in the electrostatic energy of the 
junctions. According to Josephson's law, in a non-zero voltage state, the voltage 
across the junction can be expressed as 

V. = A,.. (4) 

Each junction has a capacitance designated as Cj. Provided, the electrostatic energy 
of the i-th junction becomes 

(5) 

, and the total electrostatic energy is simply the sum, 

6 

W)=Y.'^^ii^) (6) 

i=l 

. In a straightforward manner, we have the Lagrangian: 

C{i,i) = T(f)-W(7). (7) 

Here, we apply the constraints Eq.Q by letting 

71 + 72-73 = 0, (8a) 
-72 + 74 + 75 - 76 = 0, (8b) 
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and the Lagrangian then becomes 

C = *Cf - U (7') (9) 

where, 7' = (72,73,74,75), and C is the coefficient matrix of the bihnear form for 
the charging energy. Therefore the conjugate momentum of 7 becomes 

P=^ = Cf, (10) 

where C is a tensor of rank two. By applying a Legendre transformation, we obtain 
the Hamiltonian of the system, 

n = p'f - C 

= Ip'C-'p +U{f). (11) 

Notice that here, C is a tensor relevant to the inverse mass of a certain particle" . 
From now on the system under consideration will be, given the theoretical analogy, 
referred to as the ^^particle" moving within a potential U (7') 

For a quantum mechanical treatment, one now introduces the fundamental re- 
placement, 

p = -ih^. (12) 

After this replacement the quantum mechanical Hamiltonian can be written as, 

n = Ip'c-'p 

5 

+Ei{l- cos(-72 + 73 + 27r/)} + (1 - cos 7,) 

j=2 

+E6{l-cos(-72 + 74 + 75 - 27r/)} (13) 
= r(72,---,75,/)+W(72,---,75,/) (14) 

To determine the motion of the ^^particle" , we study the reduced potential term 
W(72, ■ ■ ■ ,74, /). In Fig. El we show ^(72, ■ ■ ■ , 75, /) at / = 0.5, projected on to 
a two dimensional space. It is merely a guide to the eye, but clearly shows four 
metastable states, where the ^^partide" will be able to rest. We used a simple 
conjugate gradient method to determine the local minimums of U. The energy level 
of each metastable state can be seen in Fig. |H1 One can see that, even after a slight 
shift in /, as long as 0.5 — /c < / < 0.5 + /c, the system retains its four metastable 
states. Calculation shows that — 0.05. Each of the four states differ in current 
direction, and the behavior is symbolically expressed in Fig. ^ and the numerical 
data is hsted in Table. HI 

Flux qubits represent each of the two states, |0) and |1) by current direction. 
|0) corresponding to clockwise current in the loop and |1), counter-clockwise. The 
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FIG. 2: The reduced potential ^(72, • • • , 75, / = 0.5) projected on to a two dimensional 
space, in which the four local minimums (light-colored sections) of lA exists. Each local 
minimum corresponds to a distinct macroscopic state of the system, from the upper most 
in clockwise order: |0)|0), |1)|0), |1)|1), |0)|1). 




frustration 

FIG. 3: The energy levels of metastable states versus the frustration index /. There exits 
four metastable states when 0.5 — /c < / < 0.5 + /c, /c — 0.05. All Josephson junction 
energies are set equal to Ej. 

left and right loops will suffice to represent two distinct states, and we can see the 
system is capable of providing four quantum states, |0)|0), |0)|1), |1)|0), which 
are the basis of a two-bit quantum computer. 

The important aspect of the circuit is that it contains many parameters that 
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FIG. 4: Symbolic representation of the distinct states. The arrows indicate the direction 
of the tunneUng current. 



TABLE I: Values of the phase difference at each junction for each state, in radians. The 
symbol within the parentheses represent the tunnel current direction (through Eq.l^Q)), in 
accordance with Fig. ^ 



\m \m \m 

Junction 1 -1.52 (^) -0.57 (^) 0.57 (^) 1.51 (^) 
Junction 2 0.0 (-) 4.28 (i) 2.01 (j) O.O(-) 

Junction 3 1.51(^) 0.57 (^) -0.57 (^) -1.52 (^) 

Junction 4 -1.08 (^) 0.38 (^) -0.38 (^) 1.08 (^) 

Junction 5 -1.08 (j) 0.38 (j) -0.38 (j) 1.08 (t) 

Junction 6 1.08 (^) -0.38 (^) 0.38 (^) -1.08 (^) 



can be selected by the operator. Changes in junction energies Ei, and the junction 
capacitance Ci, varies the energy level of the metastable states. Hence fc varies 
along, and for some parameter settings, the system loses its four state configuration. 
But by adjusting the EiS and the CjS, the operator is capable of preparing optimum 
configuration, in accordance with the experimental factors. 



C. Manipulation of qubits 

The dynamical control of the qubits' state, is essential to effective quantum com- 
putation. We now introduce a scheme for manipulating the qubits in a simple yet 
efficient manner. 

We set / slightly away from 0.5, yet within the four-state regime. By this, we see 
from Fig. El that the four states each correspond to local potential minimum with 
different energies. We approximate the bottom of each local minimum of U with 
a multi-dimensional harmonic oscillator. Hence the ^^partide" wave function 
will practically be a Gaussian wave packet standing at the corresponding minimum. 
Here and on, notation may be used for and eij for its ground state energy. 

We induce transitions between states by applying a time-dependent perturbation. 
Physically, a time dependent magnetic field, oscillating at certain frequencies, ap- 
plied perpendicular to the circuit will result as a perturbation fulfilling our need. 
We assume that the particle state ket {ip) will always stay within the Hilbert space 
spanned by \i,j) = 1,2), and with this as base kets |21| , the Hamiltonian 
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plus a time dependent perturbation can be represented as, 
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^fcz represents the tunneling probability for transition 



cos {ut + . (15) 



\k,l). We have 



used, pfrom the Hermiticity of the Hamiltonian, A^^- 
introduced an approximation that, tunneling: |0,0) 



Notice that we have 
1,1) and |1,0) < — ^ |0, 1) 
are disallowed. Though this may lack rigor, by looking at Fig. |2l it should be clearly 
acceptable. The perturbation terms. 




ext 



+ 2tt{-E, sin + 73' + 2vr/) 

+ E6sin(-7^^+7r+7^^-27r/)}5/, 



(16) 



where, is the self- inductance of the k-th qubit, 50exi the amplitude of the alter- 
nating magnetic field, 6f = 6(f)ext/^o, and and 7^"' are the magnetic flux piercing 
the fc-th qubit and the phase value, respectively, while the system is in state 
Vj^ resembles the magnetic response of the system and V/j is the response of the 
Josephson junctions to the magnetic perturbation. 

Under the approximation that tunneling probability is small, we expand the 
energy eigen kets in powers of the As. The zeroth order term of the Hamiltonian 
(fl^j) contains no off diagonal elements, and therefore no tunneling occurs. By taking 
in up to first order terms, the eigen energy functions, up to a normalization factor, 
become 
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Off diagonal elements of the harmonic perturbation V cos {ut + (f) appears and in- 
terstate transitions are induced. For example. 
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FIG. 5: A Schematic diagram showing the transition process between two wave packets. 



Eq. (jl8b|l is a direct consequence of the above mentioned approximation. From 
Fermi's golden rule and time- dependent perturbation theory , transitions occur only 
between states whose energies differ from each other by fiuj. The process is described 
in Fig. El Let (t)) be the particle wave function at time t, and Cij(t) = {ip (t) 
Then, the probability amplitude evolves as follows, 

Cij [t) = Cij (0) cos I 1 j + iCki (0) e sm I 1 j , (19a) 

Cki{t) = ch(0)cos( ' — '-^t j + iCij (0) e ^'^sml ^ — '-^t j . (19b) 

Note that ip is the initial phase of the harmonic perturbation. 



1. Single bit operations 



As referred to in the introduction, single bit operations are a necessity for quan- 
tum computation. We utilize the above Eq. (fTI?|) to perform the needed operations. 
We regard the two states involved as an upstate (|t)) and a downstate(||)), and then 
we will be able to capitalize on our knowledge of the algebra of the dynamics spin 1 
systems. All unitary transformations within this two dimensinal Hilbert space can 
be decomposed into a sequence two operations, Rx (9) and Ry {6) , which are defined 
as 

(cos ~ — % sin — \ 
(20a) 
—I sm 2 cos 2 / 

/ COS I - sin § \ 
Ry (9) ^ e'^'^^y/' ={ _ I (20b) 



sm 2 cos 2 



where ax and ay are Pauli matrices. Their notation comes from the correspondence 
between them and rotation operations within a three-dimensional Euclidean space. 
Provided, physical implementations for Rx [6) and Ry (6) (referred to as rotations 
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from here and below) is sufficient for realization of all single bit operations. With 
our four state system, how do we obtain this? The solution is quite simple. 

The rotation on the ffist bit can be decomposed into two pulses. A rotation within 
the {|0, 0), |1, 0)} space and a rotation within {|0, 1), |1, 1)} will rotate the ffist bit 

successfully. From Eg. (1191) . a pulse with frequency uj = — ^ and (p = — vr with 

h 

h 

a duration time r = — j — j -6, accomplishes (6) in the ffist space, and 

2(V'o,o|V|?/'i,o) 

another pulse with frequency u = for time r = — ; — ; -6 with the 

^ 2(V^i,i|V|^o,i) 
same cp value will rotate the state ket in the remaining space. Ry (6) operations 

follows the same rule: choose the characteristic frequency, set = 7r/2, and apply 

a pulse with the appropriate area. 



Ua 



(21) 



2. Multi hit operation 

As stated in the introduction, given a complete set of unitary transformations 
for single bit operations, a CNOT implementation will, in general, be the last entry 
in the universal set. The controUed-NOT can be expressed as, in the basis {|0,0), 
|0,1), |1,0), |1,1)}, 

/I \ 
10 
1 

Vo 1 0/ 

The first bit is considered the control bit and the second the target. The greatest 
advantage of this model is that its implementation for the cnot gate is ex- 
tremely simple. A pulse with frequency uj = — — ip = it and a duration of 

n 

ttH 

T = — r- — nrrr- — r will give the desired result. 

2(^i,o|V|^i,i) 

Sections llI C ll and lll C 2l show that, by using magnetic pulses, we have a universal 
set of quantum gates for our two-bit model. Another interesting characteristic of 

our system is that, the availability of Bell states (e.g. "2 ^ ' — trivial. We 

expect this to become a trait of our system when considering quantum teleportation. 



III. DISCUSSION 



Here, we will discuss the essential aspects for our model to offer effective compu- 
tation. 
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A. Initialization 

An initial state of our model will typically be |0,0) or By setting the 

frustration index, away from the operational point (/ ~ 0.5), and letting it settle 
into its ground state. Fig. 01 shows that in the regime / ~ 0.2 or / ~ 0.8 the system 
has a definite stable state: |0,0) and respectively. From Sec. II Al we see that 
our system fulfills the initialization condition. 



B. Time scale 

The most intimidating enemy of quantum computers is always decoherence. The 
most subtle noise can ruin the whole attempt. There are two important character- 
istic times: Top and (see Sec. II A|) . We will give an order estimation of our model. 
Let the two loops have diameter of 200nm, and the Josephson junctions have junction 
areas of 200nm^ by 400nm^ hence junction energy Ej ~ 200GHz. For the circuit to 
operate in the flux regime rather than the charge regime, Ej ^ Ec = e^/C, which is 
reachable in experiment. Let Ec — Ej/lOO, the plasma frequency (eigen frequency 
of the potential bottom) Up ~ lOOGHz, self-inductance of the circuit L ~ 5pH and 
the tunnel current circulating the loop ~ lyuA. The energy difference of the states 
would be ~ 25GHz. Assume that the As in Eq.dTS} is around 0.1 ~ IGHz . 
By looking at Eq. (fTB|) . and taking ~ (po/2 ^ 10~^^Tm^, and applying a dynamic 
magnetic field with an amplitude of 1 Gauss, elements become ~ 10~^eV, and 
Vfj ~ 10~^eV. This order estimation shows that the magnetic response Vj^^ is the 
significant factor for the perturbation. Hence for a cnot sequence. Top ~ 10ns. 
Rotation operations are of about the same or less by an order. To increase the de- 
coherence time of Josephson junction circuits, many attempts have been made (l6| . 
and coherence times of up to 10/is have been observed for single junction Joseph- 
son qubits This would give us a rough estimate for our quality factor to 
be ~ 10^. However, we are not so optimistic since our system includes far more 
junctions, hence the circuit is more complicated. We therefore expect more noise 
or decoherence. Although this does not qualify the quality factor condition (see 
Sec. II A|l . we would not be overly surprised if experimental development were to 
increase the decoherence time by an order or two. 



C. Comparison with other coupling methods 



Ideas for couphng multiple qubits have been raised: placing an auxiliary super- 
conducting loop above the circuit and utilizing the mutual inductance between the 
qubit and the loop [3]. Also, directly exploiting the mutual inductance of two indi- 
vidual qubits is a possibihty. However our model has the advantage that because 
we are capable of treating the system as one element and isolating it from the envi- 
ronment all together, we expect less noise to be trapped compared to other setups. 
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There is less chance to couple to the environment and hence less noise. 
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